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p^ I SUMMARY: The Magellanic Stream (MS) is the nearest example of a gaseous 

ly^ ^ trail formed by interacting galaxies. While the substantial gas masses in these kinds 

C^ , of circumgalactic structures are postulated to represent important sources of fuel 

for future star formation, the mechanisms whereby this material might be accreted 
back into galaxies remain unclear. Recent neutral hydrogen (HI) observations have 
demonstrated that the northern portion of the MS, which probably has been inter- 
^ ' acting with the Milky Way's hot gaseous halo for close to 1000 Myr, has a larger 

VQ ' spatial extent than previously recognized, while also containing significant amounts 

' of small-scale structure. After a brief consideration of the large-scale kinematics of 

the MS as traced by the recently-discovered extension of the MS, we explore the 
, aging process of the MS gas through the operation of various hydrodynamic insta- 

\I ^ bilities and interstellar turbulence. This in turn leads to consideration of processes 

whereby MS material survives as cool gas, and yet also evidently fails to form stars. 
Parallels between the MS and extragalactic tidal features are briefly discussed with 
an emphasis on steps toward establishing what the MS reveals about the critical 
role of local processes in determining the evolution of these kinds of systems. 
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1. INTRODUCTION factor in determining the fate of the accreted mat- 
ter; e.g. does infalhng gas act as clouds on approxi- 
mately ballistic orbits or does it become part of a dif- 
, -, . , . , . r 1 I- , • fuse medium? Recent studies suggest a multi-phase 
Numerical smrulations of galaxy formation ^^^^^^ ^f the accreted material: while some inflow- 
show that galaxies can grow by accreting gas from j^g gas is shock-heated to near virial temperatures, 
cosmic filaments and satellite galaxies (e.g. Keres et a significant amount is accreted at lower tempera- 
al. 2005, Dekel et al. 2009). Even at the present tures of T < few xlO^ K where cooling times are 
time, z=0, accretion processes are likely to be ongo- short. Brooks et al. (2009) show that for a Milky 
ing and help galaxies to sustain star formation (San- ^^y (MW) type galaxy, about 60-70% of accreted 
cisi et al. 2008, Keres and Hernquist 2009). In any ^^^ jg shocked to T ~ 10^ K temperature, while both 
case, the response of the gas and any^ associated ma- ^old accretion (from cosmic filaments) and accretion 
terial on approaching a giant galaxy is an important r ■ j i u-j. /i t ^ ; „ 

^^ 6 6 6 J j:^ ^Qj. previous mergers and sateUites (clumpy compo- 
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nent) contribute each about 30-40%. While contain- 
ing a smaller amount of gas in general, the unshocked 
and 'clumpy' components play a very important role 
for building up the disk: cold gas is delivered close 
to the disk and goes on to form stars faster than the 
shocked gas, which must cool before supporting star 
formation. 

Gas stripped from satellites and by interact- 
ing galaxies are both important sources of galactic 
gas infall and provide the most direct possibilities 
for observing accretion processes. As emphasized by 
Keres and Hernquist (2009), the multiphase nature 
of galactic halos also enters into this issue by modify- 
ing gas stripping processes (see also Silk et al. 1987, 
GaUagher and Smith 2005, and Tiillmann et al. 2006 
for additional perspectives on these issues). Under- 
standing the astrophysics of satellite gas accretion 
events thus requires knowledge of the operation of 
gas stripping processes to determine the rate and na- 
ture of gas that is injected into the surroundings as 
well as models to assess the fate of the stripped gas. 

We are fortunate that the Magellanic Stream 
(MS) offers a nearby example of a gaseous remnant 
from interactions between the Magellanic Clouds 
(MCs) and the Milky Way. This feature, which ex- 
tends in an arc nearly half way across the sky, offers 
a unique, close-by laboratory to study physical pro- 
cesses of cosmological importance in the MW halo. 

In this paper, we summarize the latest obser- 
vational results and outstanding puzzles concerning 



the evolution of the MS gas. We start with recent 
observations in Section 2. In Section 3 we investigate 
the large-scale kinematics of the MS, while in Section 
4 we extensively focus on physical properties of the 
small-scale structure of the MS. We contrast the MS 
to other tidal tails in Section 5, and conclude with a 
future outlook in Section 6. 



RECENT STUDIES OF THE 
MAGELLANIC STREAM 



The MS is a huge (> 100 degree long) star- 
less neutral hydrogen (HI) structure trailing behind 
the MCs. It is a remnant of the wild past interac- 
tion of the MW with the MCs (the Large Magel- 
lanic Cloud, LMC, and the Small Magellanic Cloud, 
SMC), and of the MCs with each other. Theories 
have swung back and forth about the origin of the 
MS, with the relative importance of tidal (Gardiner 
and Noguchi 1996, Connors et al. 2006) versus ram 
pressure stripping (Moore and Davis 1994, Mastropi- 
etro et al. 2005) being still under debate. Espe- 
cially challenging in recent years have been the lat- 
est proper motion measurements (Kallivayalil et al. 
2006, Piatek et al. 2008) which increased the 3-d 
velocity for both Clouds, from 220 to 370 km s^^for 
the SMC, and from 293 to 350 km s^^for the LMC. 




Fig. 1. A collage of HI observations of the MS. The left-hand side image shows the HI column density 
distribution from Putman et al. (2003) and was produced from observations obtained with the Parkes radio 
telescope (angular resolution of 15' ). The right-hand side image shows the HI velocity field at the tip of the 
MS; this image is from Stanimirovic et al. (2008) and was obtained using Arecibo observations (angular 
resolution of 3.5' ). The color bar corresponds to the Arecibo image and shows a velocity range from —320 
km s~^ (blue) to —440 km s^^(red). 
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After incorporating these new measurements 
into the orbit calculation, Besla et al. (2007) sug- 
gested, contrary to all previous studies, that the MCs 
are not likely to be bound to the MW and may be ap- 
proaching the MW for the very first time. This leaves 
essentially very little time for the SMC-LMC-MW in- 
teractions and the formation of the MS either though 
tidal or ram pressure forces. In the mean time, sev- 
eral observational studies have suggested that the 
MW's rotational velocity is significantly higher (254 
km s~^) than the lAU standard of 220 km s^^(Reid 
et al. 2009). This implies that the MW itself is more 
massive, with a virial mass of 1.5 x 10^^ Mq. Shattow 
and Loeb (2009) showed that with a more massive 
MW binary orbits of the MCs are again possible and 
the MCs may be gravitationally bound by the MW. 

While apparently starless, the MS is the host 
of frequently detected Ha emission. The origin of 
this Ha emission has been mysterious as the expected 
Ha signal excited by the cosmic and Galactic UV 
background is significantly lower than what is ob- 
served (Bland-Hawthorn and Maloney 1999). Addi- 
tional sources of ionization have been invoked, in- 
cluding shocks, magnetic fields, and/or interactions 
between the MS clouds and the hot halo gas. 

Many past and recent low-resolution HI stud- 
ies have provided illuminating insights into the struc- 
ture of the MS. A strong velocity gradient was ob- 
served from about +400 km s^^ at the location of 

the LMC (Dec 68 dcg), to -400 km s^^ at 

the tip of the MS (Dec -- 20°), the farthest away^ 
from the MCs. Observations with the Parkes ra- ^ 
dio telescope at angular resolution of 15' (Staveley- E 
Smith et al. 1998, Putman et al. 2003, Briins et al. ^ 
2005), showed interesting large-scale HI structure in .-^ 
the form of two 100-degree long interwoven filaments ° 
(see Fig. 1, left). It was thought that the filaments > 
become overwhelmed by the MW halo around Dec 2 
~ 0°, ending up in a chaotic network of small fila- § 
ments and clumps. The total HI mass of the MS is " 
2 X 10^ M0 (Putman et al. 2003). For comparison, 
the HI mass of the SMC is 4.2 x 10* Mq (Stanimirovic 
et al. 1999). 

Most recent observations highlight two impor- 
tant observational phenomena: the MS is signifi- 
cantly more extended than previously thought, and 
has a significant abundance of small-scale structure. 
Braun and Thilker (2004) suggested that the diffuse 
northern portion of the MS extends up to Dec 40°, 
while the latest high-resolution observations by the 
GALFA-HI survey (Stanimirovic et al. 2006) showed 
a highly organized structure. Instead of a chaotic 
HI distribution, Stanimirovic et al. (2008) found 
several filamentary structures extending up to Dec 
~ 30°to the north, and reaching a heliocentric veloc- 
ity of —420 km s~^. These filaments have a great 
deal of small-scale structure, mainly in the form of 
discrete HI clouds, and have distinct HI morphologies 
and velocity gradients (Fig. 1, right). Very recently, 
Nidever et al. (2010) combined all available HI data 
sets with some new Green Bank Telescope observa- 



tions and showed that there is indeed a continuous 
extension of the northern MS from the areas cov- 
ered by Parkes and GALFA-HI surveys all the way 
to Dec 40° as suggested by Braun and Thilker (2004). 

In summary, the present-day MS is at least 
40% longer, and about 10% more massive, relative 
to the MS wc knew about a few years ago. 



WHAT PROCESSES SHAPE 
THE LARGE-SCALE HI 
STRUCTURE OF THE MS? 



Attempts to reproduce the observed HI mor- 
phology and velocity gradient along the MS have had 
a varying level of success. For pros and cons of var- 
ious models please see Connors et al. (2006) and 
Besla et al. (2008). The model-predicted HI ve- 
locity gradients especially differ at the northern tip 
of the MS, making a comparison with observations 
the easiest. For example, the Connors et al. (2006) 
tidal model predicts a heliocentric velocity of ~ —400 
km s~^at the extreme north tip of the MS, while the 
latest gravity -I- ram-pressurc model by Mastropietro 

et al. (2005) predicts Vlsr 250 km s'^ In Fig. 

2 we show these two predictions with a solid and a 
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Fig. 2. The local standard of rest (LSR) cloud ve- 
locity as a function of Magellanic longitude. Crosses 
are from Putman et al. (2003) and triangles are 
from Stanimirovic et al. (2008). Three lines repre- 
sent the predicted relationship from different models: 
Mastropietro et al. (2005, dashed line), Connors et 
al. (2006, dot-dashed line), and Besla et al. (2008, 
solid line). 



dashed line, respectively. The data points repre- 
sent MS clouds cataloged by Putman et al. (2003) 
(crosses) and Stanimirovic et al. (2008) (triangles). 
The observed velocity of ^ —430 km s^^by Stan- 
imirovic et al. (2008) is clearly in a reasonable agree- 
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ment with the purely tidal predictions. However, this 
is lower than the prediction by Besla et al. (2008) of 
< —500 km s~^at the extreme MS tip (dot-dashed 
line in Fig. 2) obtained when the latest proper mo- 
tions of the MCs are taken into consideration in a 
cosmologically-motivated gravitational model. 

This comparison clearly shows that gravity 
plays a major role in the large-scale structuring of 
the MS kinematics. While the updated mass of the 
MW (Shattow and Loeb 2009) is higher than what 
has been considered in Besla et al. (2008) and may 
affect the MS velocity, secondary effects like the gas 
drag due to the interaction of the MS with the halo 
may also be required to slow down the MS and ad- 
dress details of its spatial orientation and morphol- 
ogy- 



dynamics (SPH ) can suppress development of hydro- 
dynamic instabilities due to smoothing (Agertz et al. 
2007), while N-body simulations ignore gas processes 
altogether. Grid-based modeling can capture these 
small-scale instabilities better, exemplified by simu- 
lations exploring mechanisms for excess Ha emission 
in the MS (Bland- Hawthorn et al. 2007), galaxy 
replenishment (Bland-Hawthorn 2009, Heitsch and 
Putman 2009), and high velocity clouds (HVCs) in 
the halo (Quilis and Moore 2001). It is therefore es- 
sential to observationally constrain the effectiveness 
of various hydrodynamical instabilities for the for- 
mation and evolution of small-scale structure. 



4.1. Cloud angular size distribution 



4. WHAT PROCESSES SHAPE 
THE SMALL-SCALE HI 
STRUCTURE OF THE MS? 



Recent Arecibo high-resolution HI observa- 
tions reveal a wealth of substructure in the MS gas, 
down to an angular size of ~ 3 arcmin. Many HI 
clouds have multi-phase medium, while signatures 
for interaction of gas structures with each other and 
with the background medium are also common (Put- 
man et al. 2003, Bland-Hawthorn et al. 2007). 
The complex morphology of the HI gas, and its ob- 
served physical properties, indicate that processes 
are clearly at work on these small scales that could 
affect star formation, the transfer of gas to the halo, 
and also may provide additional drag affecting MS 
global dynamics. Although these processes play a 
crucial role for gas evolution of the MS (Murray et 
al. 1993, Bland-Hawthorn et al. 2007, Heitsch and 
Putman 2009), it is still not clear exactly how they 
operate, and on what timescales. 

For example, theoretical arguments as well 
as simulations (Mori and Burkert 2001, Quilis and 
Moore 2001, Heitsch and Putman 2009) suggest 
that instabilities act on rather short timescales (~ 
100 Myrs). This implies that the MS is being con- 
tinuously replenished, and that the stripped gas 
may eventually constitute a substantial source for 
the MW's star formation in the form of infalling 
warm ionized gas (so called "warm drizzle", Bland- 
Hawthorn et al. 2007). It is puzzling, however, that 
high-resolution observations of MS HI clouds show 
characteristics that indicate stability and longevity 
rather than rapid destruction; they are often com- 
pact, some have a multi-phase medium, and there 
are regions dense enough for the onset of molecule 
formation (which we describe and discuss below). 

It is important to emphasize that the global 
numerical simulations rarely have resolution neces- 
sary to resolve physical processes on small scales. For 
example, simulations using smoothed particle hydro- 



Stanimirovic et al. (2008) produced a cata- 
log of HI clouds and their basic observed properties. 
The cloud angular size distribution, shown in Fig. 
3, peaks at about 10', while the HI column density 
peaks at about 10^^ cm~^. If at a distance of 60 kpc 
(this is the distance of the SMC), then typical clouds 

have a radius of ~ 100 pc and a HI mass of ~ 10"^ 

M0. 
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Fig. 3. Histogram of cloud angular size (in 
arcmins) for the MS cloud population from Stan- 
imirovic et al. (2008). The dashed line shows the an- 
gular resolution limit of the Arecibo radio telescope. 



To investigate the importance of various hy- 
drodynamical effects on evolution of the MS gas we 
estimate their approximate timescales. 

(i) Thermal instability (TI) develops due to 
gas cooling and would fragment a warm stream of 
gas left behind the MCs. Assuming that the origi- 
nal MS had properties similar to those found in the 
outskirts of the SMC, the TI fragmentation will oc- 
cur on timescales < 100 Myrs (for details please see 
Stanimirovic et al. 2008). 
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(ii) Kclvin-Helmholtz instability (KHI) occurs 
at the interface between the moving warm stream 
and the hot ambient medium and provides a con- 
tinuous stripping mechanism for the MS. The KHI 
timescale depends on the properties of the halo gas 
as well (temperature, density), which are not well 
constrained observationally. However, its typical 
timescale is in the range of a few hundreds to a few 
thousands of Myrs. 

(iii) The small fragments made by TI and/or 
KHI are subject to the heat transfer from the much 
warmer ambient medium. If undergoing classical 
evaporation by heat conduction (McKce and Cowie 
1977), the HI clouds evaporate on a long timescale, 
> 1 Gyr. In the case of turbulent mixing layers the 
evaporation timescale would decrease, while a mag- 
netic field would tend to make clouds longer lived. 

To conclude, TI and KHI must have had im- 
portant effects on the shaping of the small-scale HI 
structure over the MS lifetime (in most theoretical 
frameworks at least 1 Gyr). While undergoing evap- 
oration, the HI clouds can survive for long times, and 
therefore it may not be surprising to observe such 
clumpy morphology at the MS tip. 

If we assume that TI is the dominant shaping 
agent, then we can predict a typical size of thermal 
fragments Acooi: 






(1) 



where, A is the cooling function, k is the Boltzmann 
constant, Cg is the sound speed, and T„ and n^ are 
the temperature and volume density of the warm 
neutral medium (WNM). For the density and tem- 
perature conditions characteristic of the SMC out- 
skirts (Tw = 8000 K and nw = 5 X lO'^ cm^^)^ "^yp_ 
ical" thermal fragments should be about Acooi ^ 200 
pc. A comparison with the peak of the cloud angular 
size distribution, which corresponds to 10', suggests 
that the MS tip is at a distance of ~ 70 kpc. While 
this simple, back-of-the-envelope calculation is only 
demonstrative, it is interesting that our distance es- 
timate agrees well with the predictions from tidal 
models (e.g. Connors et al. 2006). Even more im- 
pressively, our distance estimate is in agreement with 
the recent estimate of 75 kpc based on a model by Jin 
and Lyndcn-BcU (2008). This model assumed that 
energy and angular momentum are conserved along 
the MS, and that the MS is trailing on a planar orbit 
around the Galactic center. 



4.2. Multi-phase medium 



Another interesting phenomenon is the multi- 
phase HI structure of the MS. About 15% of clouds 
in the sample of Stanimirovic et al. (2008) have 
velocity profiles whose fits require two temperature 
components. This suggests the existence of a multi- 
phase medium at a significant distance from the MW 



plane. We find evidence for warm gas, with a veloc- 
ity FWHM of about 25 km s~^, and a cooler com- 
ponent, with a FWHM generally in the range 3-15 
km s^^. Similarly, Karberla and Haud (2006) inves- 
tigated velocity profiles along the MS based on the 
Leiden/ Argentine/Bonn data (Kalberla et al. 2005). 
They found that 27% of MS profiles at positive LSR 
velocities (close to MCs), and 12% of profiles at neg- 
ative LSR velocities, require two temperature com- 
ponents. 

In addition, Matthews et al. (2009) detected 
the first HI absorption lines against radio background 
sources in the direction of the MS close to the MCs. 
The two detected absorption features have a veloc- 
ity FWHM of 4.2 and 5.0 km s'^ respectively. The 
spin temperature of the absorbing clouds of 80 and 
70 K was derived, resulting in the HI column den- 
sity of ~ 2 X 10^" cm^^. The only direct detection 
of H2 in absorption is by Richter et al. (2001), who 
used Far Ultraviolet Spectroscopic Explorer (FUSE), 
and found an excitation temperature of ^ 140 K and 
7V(H2) = 3 X 10^6 cm^2_ 

Clearly, the MS contains a multi-phase 
medium. This is exciting as some cold cores are 
reaching temperature and HI column densities usu- 
ally required for molecule (CO) formation, providing 
potential for future star formation. One question 
that remains is whether these cold cores were formed 
in the MS, or have been stripped from the MCs. 
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Fig. 4. Histograms of the peak HI column density 
(in 10^* cmr"^ ) for HI clouds with a single velocity 
component (solid line) and HI clouds with the multi- 
phase structure (dashed line). The 3-a sensitivity 
limit of the GALFA-HI survey used for this study is 
-3 X IQis cm-"^. 



In Fig. 4 we show the HI column density of 
single-phase and multi-phase clouds at the tip of the 
MS. While single-phase clouds peak at A^(HI) ~ lO^^ 
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cm ^, the multi-phase clouds have higher column 
densities N(m) - 1.5 - 4 x lO^^ cm'^. The two 
distinctly different hystograms suggest intrinsic dif- 
ference between single- and multi-phase clouds, with 
the latter one being found in better shielded, more 
condensed regions. 

The level of turbulent motions of colder cores 
with respect to the warmer envelopes can be gauged 
by calculating the sonic Mach number of cold cores: 
M = \Vc- Kvl/FWHMw, with K and Kv being ve- 
locity centroids of the cold and warm cloud compo- 
nents, and FWHMw being the velocity FWHM of 
the warm component. The histogram of M values is 
shown in Fig. 5; ^ 90% of data points are within 
\M\ < 1.5 suggesting subsonic or mildly supersonic 
motions. For comparison, Heiles and Troland (2003) 
found supersonic internal motions for the MW cold 
neutral medium (CNM) clouds with M ^ 3, while 
Kalberla and Hand (2006) found that most HVCs 
have M ~ 1.5 for cold cores relative to their warm 
envelopes. 
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phase medium are clearly needed. For example, it is 
becoming well-accepted that dynamical events such 
as collisions between turbulent flows can initiate fast 
cycling of interstellar gas between multiple phases. 

While our discussion above focused only on 
the HI multi-phase structure, further insights into 
the rich multi-phase structure in the MW halo have 
been provided by UV and optical absorption stud- 
ies probing various ionization states. Detections of 
OVI absorption from gas associated with the MS by 
Sembach et al. (2003) give strong support for the 
existence of an ionized component around the MS 
with T < 10^ K. It is generally interpreted that this 
component represents an interface between the hot 
halo gas (T - 10*^ K) and the cooler MS gas. 

Studies of lower ionization states suggest the 
existence of diffuse envelopes of somewhat cooler, 
partially ionized gas that is not visible in current 
HI surveys. Specifically, Silll provides a more sensi- 
tive probe of ionized components of the infalling gas, 
likely probing different phases than OVI with T = 
-|^q4-4.5 j^ g^j^^ ^Yie corresponding TVhi ~ 10^^ cm~^, 
well below current HI survey sensitivity (Shull et al. 
2009). Such gas has been detected along a sighthne 
in the northern MS with velocities associated with 
the MS (Collins et al. 2009). Fox et al. (2010, 
in prep) have performed an analysis of many low- 
and high-ion species in UV and optical absorption, 
including Silll, against background sources in the 
northern MS. From this work, a picture emerges of a 
diffuse, multi-phase transition structure between the 
warm, mostly neutral envelope gas detected in HI 
and the hot, mostly ionized envelope gas detected in 
OVI. 



4.3. Turbulence as a generator of small-scale 
structure? 



Fig. 5. Histogram of the sonic Mach number of 
MS cold HI cloud cores relative to warm cloud en- 
velopes. 



The existence of a multi-phase medium in the 
MW halo at a distance of ~ 60 kpc is however not ex- 
pected from a theoretical point of view. Based on the 
consideration of cooling/heating processes, Wolfire 
et al. (1995) suggested that pressure-confined multi- 
phase clouds should not be found in the MW halo 
at distances larger than 25 kpc. The thermal pres- 
sure required for the co-existence of the CNM and 
WNM at a distance of 45 kpc from the MW plane 
is > 100 K cm^^, while a plane-stratified isother- 
mal MW halo with T = 10® K provides a thermal 
pressure of < 100 K cm^''. While Sternberg et al. 
(2002) showed that dark matter can provide an addi- 
tional confinement mechanism, thereby allowing the 
existence of multi-phase structure at distances < 150 
kpc, the amount of dark matter associated with the 
MS is unknown. Re-consideration of halo properties 
and/or requirements for the existence of the multi- 



Instead of being predominantly formed out of 
the smooth diffuse WNM, the clumpy and multi- 
phase small-scale structure in the MS could be a re- 
sult of turbulent inhomogeneities that originated in 
the MCs and were simply stripped during the MW- 
MCs interactions. Also, an alternative to hydrody- 
namical instabilities being the main driving source 
(as considered in previous sections), are large-scale 
shearing and tidal flows that can induce turbulence 
on large scales and provide an energy cascade and 
formation of structure on smaller scales. We briefly 
explore these possibilities in this section. 

It has been shown that the HI distribution in 
the SMC and the LMC is turbulent and can be de- 
scribed with a spatial power spectrum P{k) oc k~'^ 

(Stanimirovic et al. 1999, Elmegreen et al. 2001). 
The power-law slope of the density field is 7 = 3.4 
in the case of the SMC (Stanimirovic and Lazarian 
2001), and 7 = 3.7 for the LMC (Elmegreen et al. 
2001). One puzzling issue, however, is that these 
power spectra do not show significant changes at the 
largest sampled scales. This could be interpreted as 
evidence for interstellar turbulence being driven on 
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scales larger than the size of the SMC/LMC, > 4 — 5 
kpc. 

Recently, Burkhart et al. (2010) developed a 
new method to gauge the spatial variations of tur- 
bulence by applying the 3rd and 4th statistical mo- 
ments (or skewness and kurtosis) on the observed HI 
column density distribution. Based on MHD simula- 
tions, these high-order statistical moments are well 
correlated with the sonic Mach number. Therefore, 
by measuring skewness and kurtosis for the observed 
data, we can use the correlations derived from sim- 
ulated data sets to retrieve the spatial distribution 
of the sonic Mach number, which provides an esti- 
mate for the local level of turbulence. Burkhart et al. 
(2010) applied this method to the HI column density 
image of the SMC and found that regions with the 
highest level of turbulence are located at the bound- 
aries of the SMC bar. This suggests that large-scale 
motions between the bar and the surrounding diffuse 
HI, possibly induced by tidal interactions between 
the SMC, the LMC and the MW, or some kind of 
shearing flows, may be imprinting a strong energy 
signature on the HI gas. 

Similarly, Goldman (2000) suggested that the 
HI turbulence in the SMC was induced by large-scale 
flows from tidal interactions with the MW and the 
LMC about 2 x lO'* yrs ago. Such large-scale bulk 
flows could have generated turbulence through shear 
instabilities. If shearing flows were able to leave such 
strong imprint on the SMC, they must be also affect- 
ing the MS gas as well and may be responsible for 
the formation of the small-scale structure we observe 
in HI. 

The need for an initially clumpy MS gas has 
also been highlighted recently by Bland-Hawthorn 
et al. (2007) who proposed a shock-cascade process 
to explain the observed Ha emission along the MS. 
Two most important aspects of this study are: an 
initially clumpy distribution of the MS gas, and a 
strong interaction between the MS clouds and the 
MW halo which drives the coUisionally excited Ha 
emission. As the MS clouds upstream experience 
gas ablation by the oncoming hot Galactic wind, the 
ablated gas is slowed down and transported behind 
the clouds. The ablated gas further collides with the 
clouds downstream, resulting in shock ionization of 
HI clouds. This shock-cascade model can explain 
measured Ha intensities along the MS. The shock- 
cascade model predicts that large changes should 
take place in the HI distribution on timescales of 
100-200 Myrs. The ablation process erodes the low 
density HI gas, slowly eating into the high-density re- 
gions. As a result, after about 200 Myrs, the HI col- 
umn density distribution is highly asymmetric (see 
Fig. 6): it peaks at Af(HI) ~ 10^^ cm~^ and is miss- 
ing both low- and high-density gas relative to the ini- 
tial HI distribution. As the tip of the MS has been 
exposed to the halo the longest, the shock-cascade 
process should be clearly noticeable here. 

In Fig. 6 we compare the observed HI col- 
umn density probability density function (PDF) with 
the same quantity at two snap-shots in the Bland- 



Hawthorn et al.'s simulation: 70 and 270 Myrs af- 
ter the initial exposure of the MS to the halo wind 
(shown as dashed and dot-dashed lines in Fig. 6). 
The observed PDF was derived by taking the data 
from Stanimirovic et al. (2008), deriving the HI col- 
umn density distribution, and simply dividing this by 
5 to account for the difference in the areas probed by 
observations and the simulation (the simulated area 
is about 5 times smaller than that probed by obser- 
vations). The large difference in the simulated data 
after 200 Myrs is clearly visible, and the later distri- 
bution is missing both low- and high-density gas. 
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Fig. 6. Evolution of the HI column density PDF 
in the shock- cascade simulation by Bland- Hawthorn 
et al. (2007). The blue dot- dashed and green dashed 
lines shows times stamps in the simulation at 70 and 
270 Myrs, respectively. The solid line shows the HI 
column density PDF derived by using observations 
from Stanimirovic et al. (2008). The observed PDF 
was divided by 5 because observations sample five 
times larger projected area than the simulation. 



However, the observed PDF is not similar 
to any of the simulated PDFs. Contrary to a 
highly asymmetric simulated A^(HI) PDF, the ob- 
served PDF is highly symmetric and almost Gaus- 
sian. It clearly contains more low- and high-density 
gas than the end point of the simulation. As shown in 
Burkhart et al. (2010), subsonic turbulence produces 
Gaussian column density PDFs, while in the case of 
supersonic turbulence PDFs are highly skewed. This 
again highlights the difference between observations 
and the simulation: simulated distributions appear 
significantly more turbulent than what observations 
show. This suggests that the proposed ablation pro- 
cess is too fast and something must be slowing it 
down and helping the MS clouds survive longer. 

The structure of the boundary between clouds 
and the hot atmosphere of the MW is one factor 
which could account for the apparently slow rate of 
mass ablation in the MS. Numerical models by Vieser 
and Hensler (2007) indicate that while conductive 
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heating promotes the evaporation of a cooler, dense 
cloud moving with respect to a rarified hot medium, 
it also can significantly reduce the action of the KTI. 
Since under these circumstances KTI is likely to be 
the dominant mode of cool cloud disruption, the net 
effect of conductive heating then is to substantially 
extend cloud lifetimes. This type of mechanism also 
is consistent with the presence of gas with a range of 
ionization potentials in the MS, and merits further 
exploration. 



5. THE MS AS A VEHICLE FOR 

UNDERSTANDING CONDITIONS 
FOR STAR FORMATION IN 
TIDAL TAILS IN GENERAL 



The evolution of gas-rich tidal debris involves 
processes operating on a variety of density, tempera- 
ture and mass scales. Among these star formation is 
one interesting end point for small scale structures. 
While massive tidal tails, such as those in some of 
the Arp interacting galaxies, show clear evidence for 
active star formation (e.g. Schombert et al. 1990, 
Gallagher et al. 2001), other gas-rich tidal features, 
such as the eastern HI arm of M51 (Rots et al. 1990), 
remain as purely gaseous structures. However, the 
increased sensitivity of recent observations are re- 
vealing star formation under a wider range of condi- 
tions than previously recognized. 

On larger scales, these include the production 
of relatively large concentrations of gas that can be- 
come tidal dwarf galaxies and are capable of sup- 
porting extensive star formation (e.g. Mirabel et 
al. 1992, Bournaud and Due 2006, Weilbacher et al. 
2000, Smith et al. 2010 and references therein). A 
key point beyond this was the recognition that even 
relatively diffuse extragalactic tidal gas systems can 
contain significant amounts of molecular gas (Walter 
and Heithausen 1999, Braine et al. 2001, Taylor et 
al. 2001). Perhaps it then is to be expected that 
star formation on small scales in tidal debris also is 
showing up, especially in deep images obtained in the 
optical or in the ultraviolet with Galaxy Evolution 
Explorer (GALEX) (e.g. de Mello et al. 2008, Werk 
et al. 2008, Thilker et al. 2009). 

The MS, unlike the Magellanic bridge (Mizuno 
et al. 2006, Nishiyama et al. 2007), is not known 
to contain classical molecular clouds or candidate 
star forming regions. Thus the MS appears to be 
a large scale example of an HI stream that is sterile 
against star formation. How then does the MS differ 
from extragalactic HI features, including many tidal 
tails, that support star formation? Can studies of 
the MS provide insights into why the observed gas 
clumps evidently do not grow and become gravita- 
tionally unstable and yet also survive for long times? 
The low mean column density is likely to be one 
factor and dust content may be another. As dis- 
cussed by Maybhate et al. (2007) and Boquien et 



al. (2009), star formation usually is observed when 
N{H) > 3 X 10^" cm'^, which is larger than the 
N{HI) seen in most of the stream. 

However, since the gas in the MS can be stud- 
ied in considerable detail, it should be possible to 
go beyond this type of important but empirical bulk 
diagnosis, and carry out detailed analysis of the evo- 
lution of typical MS clumps. For example, as we 
discussed earlier, the Bland-Hawthorn et al. (2007) 
provides useful initial predictions for the evolution 
of the column density distribution within the MS. 
Comparisons between even more sophisticated mod- 
els and the new MS observations can be expected to 
yield insights not only into the astrophysics control- 
ling the MS, but more generally into the fate of HI 
injected into the vicinities of galaxies by interactions 
or other cosmologically related processes (e.g. Keres 
and Hernquist 2009). 



CONCLUSIONS, ON-GOING 
AND FUTURE STUDIES 



With the improved resolution and sensitivity 
of radio telescopes, abundant and rich small-scale 
HI structure in the MS is being revealed even in 
regions located the farthest away from the parent 
MCs and deeply embedded in the hot MW halo. The 
HI clouds often show multi-phase signatures and ap- 
pear shielded from high turbulence caused by various 
types of gas stripping by the MW halo. Occasionally, 
even cold cores with column densities that could sup- 
port molecule formation are found. 

What physical processes produce such rich 
structure in a tidal tail like the MS, how will this 
structure evolve as it interacts with the surrounding 
hot medium, and will it eventually infall to the MW 
disk? These questions have led us to explore the 
importance of various hydrodynamical instabilities 
and their effectiveness. Our analytical consideration 
of timescales, as well as recent numerical advances, 
suggest that thermal and Kelvin-Helmholtz instabil- 
ities operate on timescale much shorter than the MS 
formation time and hence must be important. Indi- 
rect studies of turbulence in the MCs suggest that 
large-scale shearing/tidal flows may be able to drive 
turbulence and cascade to smaller scales. Yet, the 
observed HI column density distribution and highly 
rich temperature structure observed in the MS over 
a range from ~ 100 to ^ 10^ K, paint a picture 
of stable, long-lived environments. One promissing 
solution could be the nature of boundary regions be- 
tween the MS clouds and the MW halo. 

To investigate the role of the MW halo in 
shaping of the small-scale MS structure, we are in the 
process of placing constraints, in the radio and opti- 
cal regimes, on characteristics of the gas in the tran- 
sition region from the MS gas to the ambient MW 
halo medium. Using the National Radio Astronomy 
Observatory Green Bank Telescope we have obtained 
the most sensitive HI emission images of portions of 



THE MAGELLANIC STREAM 



the MS to date and have begun to analyze the struc- 
ture and kinematics of the gas. Prehminary resuhs 
from this study can be found in Nigra et aL (2009). 

In the near future, a large step forward in 
understanding the MS properties will come from 
the Galactic Australian Square Kilometer Array 
Pathfinder GASKAP project. GASKAP is a study 
of the 21-cm line of HI and the 18-cm lines of OH in 
the Galactic Plane and the Magellanic Clouds and 
Stream using a new radio interferometer (ASKAP, 
Johnston et al. 2007) under development in Aus- 
tralia. ASKAP will consist of 36 12-m antennas, each 
with a focal plane phased array, and operating over 
a frequency range 700 MHz to 1.8 GHz. This new 



instrument is expected to become operational in late 
2012. GASKAP images will have an order of mag- 
nitude higher angular resolution (~ 20") and sensi- 
tivity {^ 0.04 K) relative to any previous large-scale 
survey of the MS and will constrain various eroding 
processes shaping the MS. 



Acknowledgements - It is a pleassure to thank 
Joss Bland-Hawthorn for stimulating discussions and 
sharing of simulations. We acknowledge support by 
the NSF grant AST-0908134 to SS and AST-0708967 
to JSG. SS also thanks the Research Corporation for 
Science Advancement for their support. 



REFERENCES 



Agertz, O. et al.: 2007, Mon. Not. R. Astron. Sac, 
380, 963. 

Besla, G., Kallivayalil, N., Hernquist, L., Robertson, 
B., Cox, T. J., van der Marel, R. P. and Al- 
cock, C: 2007, Astrophys. J., 668, 949. 

Besla, G., Kallivayalil, N., Hernquist,, L., vander 
Marel, R. P., Cox, T. J., Robertson, B. and 
Alcock, C: 2008, lAU Symp., 256, 119. 

Bland-Hawthorn, J. and Maloney, P. R.: 1999, As- 
tron. Soc. Pacific. Conference Series, 166, 
212. 

Bland-Hawthorn, J., Sutherland, R., Agertz, O. and 
Moore, B.: 2007, Astrophys. J., 670, L109. 

Bland-Hawthorn, J.: 2009, lAU Symp., 254, 241. 

Boquicn, M. et al.: 2009, Astron. J., 137, 4561. 

Bournaud, F. and Due, P.-A.: 2006, Astron. Astro- 
phys., 456, 481. 

Braine, J., Due, P.-A., Lisenfeld, U., Charmandaris, 
v., Vallejo, O., Leon, S. and Brinks, E.: 2001, 
Astron. Astrophys., 378, 51. 

Braun, R. and Thilker, D. A.: 2004, Astron. Astro- 
phys., 417, 421. 

Brooks, A. M., Govcrnato, F., Quinn, T., Brook, 
C. B. and Wadsley, J.: 2009, Astrophys. J., 
694, 396. 

Briins, C. et al.: 2005, Astron. Astrophys., 432, 45. 

Burkhart, B., Stanimirovic, S., Lazarian, A. and 
Kowal, G.: 2010, Astrophys. J., 708, 1204. 

Collins, J. A., Shull, J. M. and Giroux, M. L.: 2009, 
Astrophys. J., 705, 962. 

Connors, T. W., Kawata, D. and Gibson, B. K.: 
2006, Mon. Not. R. Astron. Soc, 371, 108. 



de Mello, D. F., Smith, L. J., Sabbi, E., Gallagher, 
J. S., Mountain, M. and Harbeck, D. R.: 2008, 
Astron. .J., 135, 548. 

Dekel, A. et al.: 2009, Nature, 457, 451. 

Elmegreen, B. G., Kim, S. and Staveley-Smith, L.: 
2001, Astrophys. J., 548, 749. 

Fox, A. J. et al.: 2010, in preparation. 

Gallagher, S. C, Charlton, J. C, Hunsberger, S. D., 
Zaritsky, D. and Whitmore, B. C: 2001, As- 
tron. J., 122, 163. 

Gallagher, J. S. and Smith, L. J.: 2005, Astron. Soc. 
Pacific. Conference Series, 331, 147. 

Gardiner, L. T. and Noguchi, M.: 1996, Mon. Not. 
R. Astron. Soc, 278, 191. 

Goldman, I.: 2000, Astrophys. J., 541, 701. 

Hciles, C. and Troland, T. H.: 2003, Astrophys. .J. 
Suppl. Series, 145, 329. 

Hcitsch, F. and Putman, M. E.: 2009, Astrophys. J., 
698, 1485. 

Jin, S. and Lynden-Bell, D.: 2008, Mon. Not. R. 
Astron. Soc, 383, 1686. 

Johnston, S. et al.: 2007, Publ. Astron. Soc. Aus- 
tralia, 24, 174. 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., 
Arnal, E. M., Bajaja, E., Morras, R., Poppel, 
W. G. L.: 2005, Astron. Astrophys., 440, 775. 

Kalberla, P. M. W. and Haud, U.: 2006, Astron. 
Astrophys., 455, 481. 

Kallivayalil, N., van der Marel, R. P., Alcock, C, 
Axelrod, T., Cook, K. H., Drake, A. J. and 
Geha, M.: 2006, Astrophys. J., 638, 772. 



S. STANIMIROVIC, J. S. GALLAGHER and L. NIGRA 



Keres, D., Katz, N., Weinberg, D. H. and Dave, R.: 
2005, Mon. Not. R. Astron. Soc, 363, 2. 

Keres, D. and Hcrnquist, L.: 2009, Astrophys. J., 
700, LI. 

Mastropictro, C, Moore, B., Mayer, L., Wadsley, J. 
and Stadel, J.: 2005, Mon. Not. R. Astron. 
Soc, 363, 509. 

Matthews, D., Staveley-Smith, L., Dyson, P. and 
MuUer, E.: 2009, Astrophys. J., 691, L115. 

Maybhate, A., Masiero, J., Hibbard, J. E., Charlton, 
J. C., Palma, C, Knierman, K. A. and En- 
ghsh, J.: 2007, Mon. Not. R. Astron. Soc, 
381, 59. 

McKee, C. F. and Cowie, L. L.: 1977, Astrophys. J., 
215, 213. 

Mirabel, I. F., Dottori, H. and Lutz, D.: 1992, As- 
tron. Astrophys., 256, L19. 

Mizuno, N., MuUer, E., Maeda, H., Kawamura, A., 
Minaniidani, T., Onishi, T., Mizuno, A. and 
Fukui, Y.: 2006, Astrophys. J., 643, L107. 

Moore, B. and Davis, M.: 1994, Mon. Not. R. As- 
tron. Soc, 270, 209. 

Mori, M. and Burkcrt, A.: 2001, Astron. Soc Pa- 
cific Conference Series, 222, 359. 

Murray, S. D., White, S. D. M., Blondin, J. M. and 
Lin, D. N. C: 1993, Astrophys. J., 407, 588. 

Nidever et al.: 2010, in preparation. 

Nigra, L., Stanimirovic, S., Gallagher, J. S., Ill, 
Lockman, F. J., Nidever, D. L. and Majew- 
ski, S. R.: 2009, arXiv.org, 0908.4218. 

Nishiyama, S. et al.: 2007, Astrophys. J., 658, 358. 

Piatek, S., Pryor, C. and Olszewski, E. W.: 2008, 
Astron. J., 135, 1024. 

Putman, M. E., Staveley-Smith, L., Freeman, K. C, 
Gibson, B. K. and Barnes, D. G.: 2003, As- 
trophys. J., 586, 170. 

Quilis, V. and Moore, B.: 2001, Astrophys. J., 555, 
L95. 

Reid, M. J. et al.: 2009, Astrophys. J., 700, 137. 

Richter, P., Sembach, K. R., Wakker, B. P. and Sav- 
age, B. D.: 2001, Astrophys. J., 562, L181. 

Rots, A. H., Bosma, A., van der Hulst, J. M., 
Athanassoula, E. and Crane, P. C: 1990, As- 
tron. J., 100, 387. 

Sancisi, R., Fraternali, F., Oosterloo, T. and van der 
Hulst, T.: 2008, Astron. Astrophys. Reviews, 
15, 189. 



Schombert, J. M., Wallin, J. F. and Struck-MarceU, 
C: 1990, Astron. J., 99, 497. 

Sembach, K. R. et al.: 2003, Astrophys. J. Suppl. 
Series, 146, 165. 

Shallow, G. and Loeb, A.: 2009, Mon. Not. R. 
Astron. Soc, 392, L21. 

ShuU, J. M., Jones, J. R., Danforth, C. W. and 
Collins, J. A.: 2009, Astrophys. J., 699, 754. 

Silk, J., Wyse, R. F. G. and Shields, G. A.: 1987, 
Astrophys. J., 322, L59. 

Smith, B. J., Giroux, M. L., Struck, C. and Hancock, 
M.: 2010, Astron. J., 139, 1212. 

Stanimirovic, S., Staveley-Smith, L., Dickey, J. M., 
Sault, R. J. and Snowden, S. L.: 1999, Mon. 
Not. R. Astron. Soc, 302, 417. 

Stanimirovic, S. and Lazarian, A.: 2001, Astrophys. 
J., 551, L53. 

Stanimirovic, S., Putman, M., Heiles, C, Peek, J. 
E. G., Goldsmith, P. F., Koo, B.-C, Krco, 
M., Lee, J. -J., Mock, J., Muller, E., Pandian, 
J. D., Parsons, A., Tang, Y., Werthimer, D.: 
2006, Astrophys. J., 653, 1210. 

Stanimirovic, S., Hoffman, S., Heiles, C, Douglas, 
K. A., Putman, M. and Peek, J. E. G.: 2008, 
Astrophys. J., 680, 276. 

Staveley-Smith, L., Kim, S., Putman, M. and Sta- 
nimirovic, S.: 1998, Rev. Mod. Astron., 11, 
117. 

Sternberg, A., McKee, C. F. and Wolfirc, M. G.: 
2002, Astrophys. J. Suppl. Series, 143, 419. 

Taylor, C. L., Walter, F. and Yun, M. S.: 2001, As- 
trophys. J., 562, L43. 

Thilker, D. A. et al.: 2009, Nature, 457, 990. 

Tiillmann, R., Breitschwcrdt, D., Rossa, J., Pietsch, 
W. and Dcttmar, R.-J.: 2006, Astron. Astro- 
phys, 457, 779. 

Vieser, W. and Hensler, G.: 2007, Astron. Astro- 
phys., 472, 141. 

Walter, F. and Heithausen, A.: 1999, Astrophys. J., 
519, L69. 

Weilbacher, P. M., Due, P.-A., Fritze v. Alvensleben, 
U., Martin, P. and Fricke, K. J.: 2000, Astron. 
Astrophys., 358, 819. 

Werk, J. K., Putman, M. E., Meurer, G. R., Oey, 
M. S., Ryan- Weber, E. V., Kennicutt, R. C, 
Jr. and Freeman, K. C: 2008, Astrophys. J., 
678, 888. 

Wolfire, M. G., McKee, C. F., HoUenbach, D., Tie- 
lens, A.: 1995, Astrophys. J., 453, 673. 



10 



THE MAGELLANIC STREAM 



CTPYKTyPA HA MAJIMM CKAJIAMA MAFEJIAHOBOr 
nOTOKA KAO OCHOBA 3A TAJlAKTHT-IKy EBOJiyilMjy 

S. Stanimirovic, J. S. Gallagher III and L. Nigra 



Astronomy Department, University of Wisconsin, Madison, J^15 North Charter Street, 

Madison, WI 53711, USA 
E-mail: sstanimi@astro.wisc.edu 



UpezAednu pad no nosuey 



MarejiaHOB hotok (Mil) je Haj6jiH>KH npHMcp 
racHor ipara (})opMHpaHor MeljyAejcTBOM 
rajiaKCHJa. KaKO je SHaMajna Maca raca y 
OBHM BpcTaMa OKorajiaKTHMKHx CTpyKTypa 
npcTnocTaBJteHa jia, npeACTaBJBa BasKan 
H3Bop ropHBa 3a 6y/i;yRe (|)opMHpaH)e 
BBes^a, MexaHHsaM kojh OBaj MaTepHJan 
MO>Ke BpaTHTH Hasa/i y rajiaKcnjy oc- 
Taje nejacan. He^aBHa nocMaTpaiia ney- 
TpajiHor BOAOHHKa (HI) cy noKasajia ^a 
cenepHH ^eo Mil, kojh BcponaTHO HHiep- 
aryje ca bpcjihm racoM hs xajioa Mjie- 
^Hor IlyTa 6jiH3y 1000 Myr, aaxBaxa 
BetiH npocTop Hero niTO je npeTxo/iHO 
6hjio yTBpIjCHO, Ma/ia laKolje caAP>KH 
SHa^aJHy KOJiHMHHy cipyKTypa na Majioj 



CKajiH. Ilocjie KpaTKor pasMaTpaiba 

KHHeMaTHKe Ha bcjihkhm CKajiaMa y Mil 
ycjiOBJBCHOM HC/iaBHO OTKpHBenoM npom- 

HpCHOCTH Mn, MH Hsy^aBaMO CTapOCT 

npon,eca y racy MH Kpos agjctbo pasjiH^H- 

THX XH^pOZlHHaMH^KHX HCCTaGHJIHOCTH H 

MeljysBCsaaHHx TypGyjicniiHJa. Obo nac 
BO/iH j\o pasMaTpaita npoijeca y KOJHMa 
MaTcpHJaji MH npesKMBJLaBa Kao xjiaAna 
racna ^'A3'a, ta join yBCK BBTAjxewmo ne 
ycneBa ^a 6;>o\iMTAY>a. SBes^e. Ilapajie- 
jie H3MeI)y Mil h BanrajiaKTH^KHx hjihm- 
CKHx (|)opMan,HJa cy KpaiKO AncKyTOBane ca 
HarjiamaBaiBeM KopaKa kojh ycnocTaBJtajy 
niTa Mn OTKpHBa o npaBHJiHMa BCsanHM sa 
jiOKajiHC npoH,ccc y oztpdjHBaity CBOJiyiiHJe 

CHCTCMa OBaKBC BpcTe. 
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